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Abstract 

We consider newtonian dynamics of N charged particles on the circle with nearest neigbour interaction 
with Coulomb repulsive potential r^ 1 . Also there is an external accelerating force which is nonzero only 
on a small part of the circle. We construct homogeneous solutions where the velocities of all particles are 
approximately equal and their density is approximately uniform. This gives a qualitative mathematical 
model for some features of the direct electric current (DC), in agreement with a suggestion by R. Feynman. 

1 Introduction 

The most developed part of mathematical statistical physics is the theory of equilibrium (Gibbs) states on dis- 
crete lattices. This science is based on simple axioms followed by many solved problems, including very difficult. 
However, in nonequilibrium statistical physics in continuous space there are problems not even formalized on 
the mathematical level. 

One of such examples is the direct electric current (DC). On the macro level it is described by Ohm's law, 
and, on the micro level, it is often presented as the classical system of free or weakly interacting electrons, 
each accelerated by the constant external force and impeded by external media. There is an abundance of such 
models: the first such model refers to Drude, 1900, one can find it in any textbook on condensed matter physics, 
see for example pQ. Recent papers deeply investigate possible friction mechanisms, see for example [21 El SI - 
Here, on the contrary, we ignore this friction problem but turn to another fundamental problem. 

The question arises where the accelerating force comes from, because among hundred kilometers of power 
lines the external force acts only on some meters of the wire. Here what one can read in the "Feynman lectures 
on physics" (volume 2, sec. 16-2): 

"...The force pushes the electrons along the wire. But why does this move the galvanometer, which is so far 
from the force ? Because when the electrons which feel the magnetic force try to move, they push - by electric 
repulsion - the electrons a little farther down the wire; they, in turn, repel the electrons a little farther on, 
and so on for a long distance. An amazing thing. It was so amazing to Gauss and Weber - who first built a 
galvanometer - that they tried to see how far the forces in the wire would go. They strung the wire all the way 
across the city...". 

This was written by the famous physicist. However, after that, this "amazing thing" was vastly ignored in 
the literature. Many more questions arise. For example, why DC moves slowly but its stationary regime is being 
established almost immediately. Here we are occupied with the first one (some results concerning the second 
problem see in [7]). Namely, we want to demonstrate rigorously that even on the classical (non-quantum) level 
there is a mere possibility that the stationary and space homogeneous flow of charged particles may exist as a 
result of self-organization of strongly interacting (via Coulomb repulsion) system of electrons. 

We use only classical nonrelativistic physics - Newtonian dynamics and Coulomb's law, but also the simplest 
friction mechanism, ignoring where this friction mechanism comes from. Rigorous models of strongly interacting 
electron systems, interacting also with the ionic lattice, do not exist now. 

The Model We consider N point particles i = 1,2, N initially at the points 

< xi(0) < ... < x N (0) < L 

of the interval [0,L] <E R. We assume periodic boundary conditions that is we consider the circle Sl = [0,L) of 
length L. The trajectories Xi(t) are defined by the following system of N equations 

= ~dxl + ^ ~ <^ « 
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The interaction U between the particles is 



N 



U{{ Xi }) = Y,V{xi+i-Xi) 



where of course afjv+i — x\ and 



i=l 



V(x) = V(-x) = - > 0,r 
r 



The case a > corresponds to the Coulomb repulsive potential, that we consider here. Then the repulsive force 
is 

f(r) = - 



dV ( r ) _ _-a 



dr 

It follows that the particles, during the movement, cannot change their order. gF(x) - an external accelerating 
force (assumed sufficiently smooth) with scaling parameter g > 0. The friction function a(v) is specified below. 
It defines the loss of kinetic energy via the interaction with external media. 

It is well-known that the solution of the system ([1]), for any initial conditions, exists and is unique on all 
time interval [O.oo), under sufficiently general assumptions on the functions F and a. However, to get more 
detailed information about trajectories, one needs sufficient efforts. 



2 Effective force 

Parameters and constants Throughout the paper we are dealing with macro and micro parameters and 
(absolute) constants. Absolute constants do not depend on the parameters of the model. 

All our constructions are for N sufficiently large but finite - one cannot directly perform the limit N — > oo 
because there will be different micro-scales, influencing on the macro-parameters. Macro-parameters L and 
F(x) are fixed (do not depend on N), for example we put 

C(F,m) = max \F i - m \x)\ 1 m = 0,1,2, 

X 

and put for convenience 

C(F,0) = 1 

Roughly speaking, our first approximation dynamics is 

x k (t)=x k (0) + Vt,^ft = V 
at 

where V is a macro-parameter (the approximate velocity of particles). 
Micro-parameters 

M = M^ N \g = g^\a = a^ 

depend on N (they satisfy some conditions defined below), but we will omit index (N). For example, in some 
physical situation (in SI units) approximately 

AT- 1 = 10- 10 ,M = 10- 30 ,a = 10- 28 

These numbers were some guide for us, but we could not fit them completely - g had to be assumed smaller 
than necessary. 

Static configurations For any particle configuration (a;, = x[ N \ Vi = vf 1 ^) 

< x x = x[ N) < ... < x N = a$° < L 
the effective force, acting on the particle i, is 

W{ - a(vi) 

where 

Wi = w(xi-x,Xi,Xi+x) = f{xi - Xi-x) - f{xi+x - Xi) + gF(xi) 
The following crucial result depends only on the parameter C a ^ g — a~ 1 g. 



Lemma 1 Assume C a ^ g fixed or bounded as N — > oo, then for any sufficiently large N there exists configuration 
x\ < ... < xn (assuming zero velocities) such that the effective force is the same for all i = 1, ...,N, that is there 
exists w such that 

f(xi - x^i) ~ f(x i+ i - Xi) + gF(xi) = Wi =w,i = l, N (2) 
Moreover, for this configuration the following properties hold: 

1. uniformly in % = 1, N as N — > oo 

— 

— &i rNJ ^y- 

2. for any i denote A.; = A(l + <$i), Ai = A. Then for all i 

<4a- x sioA 

where Lo is the length of the support of F(x). 

Proof. Let us call ip(x) = gF(x) — w the virtual force. Then the required configuration can be interpreted 
as a fixed point of N particle system for the external virtual force. The virtual force is potential iff 



tp(x)dx = g F(x)dx — Lw = (3) 
Js L 

Then the (virtual) potential of the virtual force is 

W(x) = - [ i/j(x)dx 
Jq 

and a required fixed point exists as a global minimum (in R ) of the potential 

JY 



U(X!,...,X N ) +J2 W ( X i) 



If such minimum is not unique, we take anyone. 

The effective force w > can be found from condition ([3]). At the same time, summing up the equations 
©, we get 

N 

»=i 

Thus, the constant effective force equals 

If 1 N 

w = g-\ F{x)dx = g-Y.F^i) (4) 

Assertion 1 of the Lemma was proved in Theorem 1 of [5], see also [6]. Now let us prove the assertion 2 of 
the lemma. Summing up the equations ^ for i = 2, k, we get 



or 



Then 



/(A(l + 4)) - /(A) = J>F(xO - w) 

i=2 
k 

(1 + 5 k )- 2 - 1 = a- 1 A 2 Y.^F{xi) - w) 

oo k 

= [1 + Qfc]~* - 1 = J2 ® k = Y.^ F (^) - «0 (5) 

m— 1 i—2 

, ,. m 1.3...(2m-l) (2m)! 1 

2 m m\ (2 m m\y y/irm 



and 

k 

W<gijr, \£b,F{ Xi )-w)\<(l + ^)g{k-l) (7) 

i=1 

It follows that for sufficiently large N 

if a~ x {l + ^-)gL A < \. The Lemma is proved. 

Remark 1 The assertion 1 of Lemma[l] says that the asymptotics of Afc does not depend on F. We say in this 
case that F is not seen on the microscale TV -1 but only on the sub -micro scale, see Jf)j/ and J^. It follows that 
the density is macro-homogeneous. However, as it is clear from {5p and {6]), on the interval where F(x) = 0, 
the distances between particles slightly increase in the clock-wise direction. 

3 Main result - macro- homogeneous dynamics 

For any sufficiently large TV we shall prove existence of the dynamics for t <= [Q,T),T — T(N) sa N, which we 
call macro- homogeneous on [0, T). That is we shall prove that the following two properties hold as N — > oo: 

1. (asymptotically homogeneous velocities) There exists constant (one more macro-parameter) V > such 
that uniformly in i = 1, 2..., N and t € [0.T) the velocities 

Vi{t) = vf\t)^V 

2. (asymptotically homogeneous density) uniformly in t € [0.T), for any interval I C Sl the number N(I,t) 
of particles in / 

Assumptions We choose the simplest friction mechanism defined by the function which is linear in the (micro) 
vicinity of the macroparameter V 

a(v) = A Q + Av,v e (V - A, V + A) 

where Ao and A > are microparameters, coordinated with macroparameter V and microparameter w so that 

a(V) = A + AV = w > (8) 

Concerning the parameters, roughly speaking, there are two assumptions: g is small enough, in particular g <§C A 
and 

N^A 2 ^ 1 < M < min(aiV~ 1 , A) 

More exactly, it can be formulated as follows. For any sufficiently large N and some sufficiently small absolute 
constant p > 

N~ 1 A 2 a - 1 p- 1 < M < pmin{aN- 1 ^^, A) (9) 

C A , g = A~ x g < pN~ 3 (fO) 
In particular, we will use below the first inequality in (0) more concretely 

MaN 



> (n^rr 1 (n) 



A 2 

The following example shows that these inequalities provide non-empty and natural domain of parameters 

g = N^ 1 , a = N~ 12 M = N^ 3 
with constants 71 > 2, 72 > 0, 73 > and such that 

71 > 72, 1 + 271 - 72 > 73 > max(l + 72, 71) 

Put also 



L 

C{F\int) 



C w>g = w^g = 7 ^ rB .^ ,C{F,int) = \ F{x)dx (12) 



Initial conditions We fix the initial configuration for the required dynamics as Mj = V and Xi are chosen as 
in Lemma [TJ that is 

x k+1 (0) - x k (0) = A k ,w = aA-^ - aA~ 2 + gF{x k {0)) 

and it is convenient to choose the coordinate system by xi(0) = 0. Thus the initial effective force is zero. 
We could choose the initial velocities as 

u*:(0) =V + u k {0) 

where u k (0) could be assumed sufficiently small, but for convenience we always assume u k (0) = 0. 
Main result 

Theorem 1 Assume @) and (9))- il0\) . Then for the chosen initial conditions (explicitely defined in the next 
section) there exists f3 > 0, not depending on N, such that the dynamics is macro-homogeneous on the time 
interval [O.T),T = T N = f3N. 

Plan of the proof Choosing initial configuration x k (0) is a delicate matter. We did it in the previous section 
so that the effective forces acting on each particle are the same at time 0. This dynamics does not satisfy the 
main equations and we derive (in section 4) the equations for the deviations y k (t) — x k (t) — x k (0) — Vt. It is 
very important that our choice of the initial conditions allows to exclude the constant component of the effective 
force in the equations for y k {t). 

We fix the basic linear part of these equations for y k {t) and solve them (section 5). It is rather straightforward 
but demands delicate estimates. The rest linear and nonlinear parts of the equations are considered as the 
perturbation and demand some iteration procedure. The section 6 is devoted to convergence of this procedure 
and to stability estimates. To prove this we introduce special Banach space where the convergence holds, that 
also demands some nontrivial estimates. There are many similarities in the estimates - we tried not to repeat 
them. Finally, in section 7, we give some remarks and perspective. 

In the rigorous proof we tried to be very accurate with micro-parameters, that is with the parameters 
depending on N. Absolute constants (we will meet finite number of them) are denoted c, Cx,c%, ... and could 
be easily explicitely written but we did not do this because of no interest. Also sometimes we denote C(F) a 
generic macro-constant depending only on F. Starting from section 5 we take L = 1. 

4 Equations for the deviations 

As the force F(x) is not translation invariant, the dynamics x k (t) = x k (0) + Vt. cannot satisfy equations (fl}. 
We introduce the deviations y k (t) and their velocities u k {t) by 

x k (t) = x fc (0) + Vt + y k (t), y k (0) = 0, 

v k (t) = ^,«fc(*) = ^jf = "*(*) - V>Vk$) = J Mt)dt, 
and rewrite the main equations (fTJ) as equations for the deviations 

M^- = /(A fc _i(0) + y k (t) - y k -i{t)) - f(A k (0) + y k+1 (t) - y k (t))+ 

+ gF(x k (0) + Vt + y k {t)) - a(V + u k (t)) (13) 

For any function h(k) on the finite cyclic group {k : k = 1, ...Nl}, with zero element Nl — 0, introduce the 
shift operator 

(Sh)(k) = h(k + l) 
and finite difference operators (discrete derivatives) 

Vfcft(fc) = Vh(k) = V + h(k) = h(k + 1) - h(k), (V"/i)(fe) = h(k) - h(k - 1), 

in particular 

Vyjt(£) = y k +i(t) - y k (t) 



Sometimes we will use Leibnitz formula for discrete derivatives (concerning calculus of finite differences see (8] 
and references to classical papers therein) 

n 

v(/ x .../ n ) - (v/i)S(/ 2 .../„) + /iV(/ 2 .../„) = ... = ]T /i-/ m -i(v/ m )S(/ m+1 .../ n ) (w) 

Note that it differs from the standard Leibnitz formula for differentiation only by shift operators. As we will 
use it only for estimates from above, which are always uniform in k, the shift operators will not play role. 
Then 

/(A fc _i(0) + y k (t) - Vk-i(t)) - /(A fc (0) + y k+ i(t) - y k {t)) = a(-V-)(A fc (0) + V + y k y 2 = 

= a(-V-)(A,T 2 -2^|W)+L 3 (15) 



L 3 = a(-V~)[(A fc (0) + V + y k )- 2 Af + 2^^] 



A 3 fe 



Formally 



L 3 = a(-V-)A- 2 £ rf m ( V "'"f ( * ) ) m , d ro = (-l) m (m + 1) (16) 



A fc 

Applying -«V~ to the first and second terms in the right-hand side of (fl~5j). and putting 

4 = (i+4r 3 -i 

we have 

- a V-A^ = -a(A^-A^), 
a(-V-)[-2A- 3 V+y fc (t)] = 2a[A- 3 V+y k (t) - A^V+^t))] = 

= 2aA- 3 [(l + 4) _3 (y fc +i - Vk) - (1 + <Jfc-i) _3 (j/fc - Mb-x)] = 2aA- 3 [y fc+1 (t) - 2y k {t) + y k _i(t)] + L 2 

i 2 = 2aA- 3 V-[4Vy fc ] 



where 



Then the main equations become 

+ A^- + [aAf - aA-^ - gF(x k (Q)) +w}- 2aA^[y k+l - 2y k + y k ^)] = 

= L 2 + L 3 + gF(x k (0) +Vt + y k (t)) - gF(x k (0)) ~ A%) 

at 

As the first square bracket (that is the "effective acceleration force") is zero, we get the final form of the 
(differential-difference) equations 

M fyk + A dy^t)_ _ 2aA - 3[yk+i _ 2yk + =Li+L2 + l 3 + gMt) (i7) 

dt' i at 

where 

Mt) = F(x k (0) + Vt) - F(x k (0)) 

L x = g[F(x k (0) + Vt + y k (t)) - F(x k (0) + Vt)] 
First of all, we shall study in detail a cut-off system 

+ - 2aA^[y k+l 2y k + „ fc _ x )] = g^ k {t) (18) 

which we call the basic linear approximation and consider Lj, i = 1, 2, 3 as perturbation terms. We take them 
into account in section 6. 



5 Linear stability 

From now on we put L = 1 , N 



— N\. Denote yk,o(t) the solution of equations ([T8]) . 



5.1 Fourier transform 

Denote 

= $(/i)(n) = ^ exp(27rm|:) (19) 

fc=i 

the Fourier transform 



iV 



$ : ; 2 ({ A : fc = l, TV}) ; 2 ({ n = , 1, ...,N- 1}) 
of the function h = hk on : fc = 1, N}. The inverse Fourier transform is 

h k =J2 ex P (-2^nA) (20) 

n 

Put 

77(71,*) = $(Mk,o(t)), = $ (<M*)) 

Multiplying (| 18|) on A- exp(27rm-^), summing in k and dividing by M, we get the decoupled equations for 
n = 0,l,...,JV-l 



d 2 ri(n,t) A dn(n,t) 4a, „ 1 . , . „ , -i ~, . 



The characteristic equation 



Q(z) =z 2 + ^z + ^A- 3 (l - cos27rnA) = 



has the roots 



If 



A , / 16Mq , 1 
*i, 2 (n) = - — (1 ± a/1 - -^(1 - cos 27m-)) (22) 



(21) 



1 -M (1 - CO82ra iv )/0 

then the roots are different and the general solution is 

r)(n,t) = r)i(n,t) + r) 2 (n,t) 

where for I = 1 , 2 

rft(n,i) = Q(7i)e Z!t + ? ^ T / e'^M^g^t^dh = 
W \ z i) Jo 

tut— 1 /•< /•< 

= Q(n)e Z!t + - T / e^-^g^nM)^! = C,(n)e"* + r,(n)fl / ^^^(n.ti)^ (23) 

2;?; + Jo Jo 

where 

n(n) = (~l)' , 1 (24) 

| V / l-^f(l-cos27rni)) 

As j/fc(0) = we have 

Introduce the Fourier transform of Uk (0) 

r?'(n,0) = ^|^= C l {n)z l {n) = C 1 {n){z l -z 2 ) 

1 = 1,2 

It follows 

\C 1 \ = \C,\ = \^-\<^\r 1 '{nM 
z\ — z-i A 

From now on we assume for simplicity that Ufe(0) = for all k, then 



5.2 Main "linear" lemmas 

If 



-(2tt) 2 > 2 



A 2 A 

(garantied by the left inequality of ([5])) then for any n ^ the roots Z\^{n) are complex conjugate. 

In the linear case we need only boundedness of C A ,g and C' a ^ g and condition (JSJ. Under these conditions we 
will prove 



Lemma 2 



and for any n^O 



|T?(0,t)| < A 2 (t + l)C A , g (C(F)+4C a , g ) (25) 



Proof. We have for n ^ 



KM)I<-CaW— (26) 
n \ a 



n 7 2M' 



As, by the left inequality in ©, is sufficiently large, and 1 — cos27rn-^ > ci(-^) 2 for some c\ > and any 
n^O, then for I = 1, 2 from (pM|) we have 

0|r,(n)| = 2C Ag ] . 1 1 < ~ ( A ' 



^1-^(1-008 2^)) ^/^l(f) 2 ci 



Then 



|^(„, t )| < 2Ca ' 9 I f'e^-^dh] = . 2Ca ' 9 jjR^r^l - e Rzit ) < 



K)A/g ^n^2„ Jo / 16Ma ^ti\2„ 

A 2 A 3 \ AT ' L l V A 2 A 3 vat J c l 



<cW— (^) 
n V a 



^(0) = -^,^(0)=0,|r lj3 (0)| = | (28) 



For n = we have by (|55J and (|2"3|) 

and we need the following 
Lemma 3 For n = we have 

sup |0(O, i) | < A 2 C{F,2)+8C(F,l)C a . g A (29) 
t 

I f 4>(Q,t 1 )dt 1 \<A 2 (2(t + l)(C a . g + C(F,2)) + 2C(F,l)) (30) 
Jo 

From (|28[) and this lemma we get 



M0,i)| = 5 |r 1 (0)||^ 1 (0)|- 1 sup |^(0, t)| <2C A ^(A 2 C(F,2)+8C(F,l)C a , g A) < 



< 2C A J^A 2 C(F, 2) + 8C(F, 1)A 2 ) 



as ^C a g < A by ©. Also 



A ^ a >9 

\ m (0,t)\=g\n(0) f j>{Q,t x )dh] ' ^' 

/n 



< 2C^ g A 2 (2(i + 1)(C . + C(F, 2)) + 2C(F, 1)) 



1*7(0, t)| < \m(0,t)\ + \ m (0,t)\<A 2 (t + l)C A JC(F)+4C a , g ) 

where 

C(F) = 6C(F, 2) + 24C(F, 1) 

Now we can prove the following 
Lemma 4 Uniformly in k 



|lfe,o(*)l < cCa, s \ —lnN + A 2 (t + l)C A , g (C(F) + 4C a . g ) 
1 a 



|V 2 y fc+ i, (t)| < cC A , g \ —InN 

a 



This is the direct calculation via Fourier transform. We use inverse Fourier transform (f20|) , thus we have 
only to sum up the terms of Lemma [2j this gives 



\!Jk 



oWi<E E fowl 



We see that the main term corresponds to the zero mode. Noting that 

*(Vifc,o(t)) = V(n, t)(e~ 2 ^ - 1), <i>(V 2 y M (i)) = r,(n, t)( e - 2 ^ - l) 2 (31) 

and 

\Vy k , (t)\ < E E h(n,t)(e- 2 ^ - 1)1 (32) 

n^O 2=1,2 

|VV,o(t)| < E E I^M)(^ 2 " f - 1) 2 | (33) 

ri^O 2=1,2 

We see that for the first and second differences the zero mode vanioshes, that gives better estimate. The Lemma 
is proved. 

5.3 Proof of Lemma [3] 

To prove (|2"5)) note that for any t there is m(l,t) such that 

(0) <X!(0) + Vt < x m(ltt)+1 {0) 

and denote 

7 = an(0) + Vt - a; m( i,t)(0) =Vt- x m{ljt) (0) 
Put m(k,t) = m(l,t) + (k - 1). Then 

a!fc(0) + Vt-a: m ( fc) t)(0) =i + m (34) 

and uniformly in k 

\m\ <8C Q:9 A 

This follows from ([5]) and the evident formulas 

fc-i fe-i fc-i 
ar fc (0) = a:i(0) + E Ai = E A (l + «*) = (fc - 1) A + A E ^ 

i—1 i—1 i—1 



m(k,t)-l 

x m (k,t) (0) = x m{l , t) (0) + (fc - 1) A + A ^2 5 i 

i—m(l,t) 

m(fc,t)-l 

IEfc(0) + 74-:E m(M) (0)=7 + A( ^ Si-J2 § i) 

i=ro(l,t) i=l 

Then we can write 

k k 

= Jj E fW (%(M) (0) + fc)(7 + %) 
fc 

where \0f-\ < 7 + |?7fc| and we used But 

l^FW(x m(M) (O) + fc )^A^oo / F( 1 )(x)dx = 

and moreover 

|l^F( 1 )(a: m(feit) (0)+^)|7=|^^F«(x m(M) (0) + ^)- f F^da^ < 

< C(F, 2)A 7 < C(F, 2)A 2 

|^^F (1) (z m (fc, t )(O)+0 fe >7fc| < C(F, l)8C Q , ff A 
fc 

Finally the bound is 

C(F,2)A 2 + C(F,l)8C Q , g A 

To prove (|30|l is more difficult: one should obtain maximal cancellation by carefully grouping the summation 
and integration terms. Denote 

b m = min(A m _i, A m ),I m = (x m (0) - -^,a: TO (0) + ^) C Si 
For any pair (k,m) and any 4 define the set T(fc, m) = T(k,m,t) C [0,4] as follows: t\ <E T(k,m) iff 

x k (0) + Vti el m ,h <t 

Note that for given k the sets T(k, m), to— 1, N, are disjoint, similarly for given m the sets T(fc, m), fc— 1, iV, 
are disjoint. The set T(k, to) consists of disjoint intervals 

Ji = Ji(k,m),i = 1,2, ...,/3 = /3(fc, to,4), 

enumerated in the order of their hitting (imagine a particle starting at x&(0) and moving with constant speed 
V), and /3 — [tV] + l or (3 = [tV] + 2. All these intervals have length |«7j| = V~ 1 b m , except Ji(fc, fc), having 
length 7}V~ 1 bk (we call them initial intervals), and, for given fc, possibly one of the others (namely Jp, we call 
them end intervals), which length can be less than V^ 1 b m . We call the intervals Ji(k,m), having length fe, 
regular, the others - non-regular. 

Lemma 5 The sets Ti(k) = [0,4] \ U m T(fc,rn),fc = 1, ...,N, and T 2 (m) = [0,4] \ U fc T(fc, to), to = 1, ...,7V, have 
measure less than 2(4 + l)C ai9 A 2 . 



In fact, from ([5]) it follows that 



S m - 5 m -! = ai(Q l m - Q l m _j) = -a- 1 A 2 (gF(x m ) - w) + ^ MQL ~ QL-i) 

1 = 1 1=2 

\S m -5 m ^\<2C a , g A 2 (35) 

Then the assertion follows from 

N 

\A m - A m _i| = A\6 m - 5m-i\ < 2C Q , ff A 3 , ^ l A »> - A ™-il < 2C Q , S A 2 

m— 1 

Lemma [S] is proved. 
We can write 

r t i N i r l N N 

/ 0(O,ti)tfti= / dt 1 -J2(F(x k (0) + Vt 1 )-F(x k (0))) = - dh(£ / F(x k (0) + Vt 1 )-J2F(x m (0))) = 

■* •'° fe=l ,/o k=l m=l 

. N N P f r 1 JV 

^BEE +/ )n^0) + W 1 )dt 1 --£^(* m (0)) (36) 



iV 

and for any m 



AT /9 

tF(x m (0)) = F(x m (0))Q2^2\Ji(k,m)\ + \T 2 {m)\) 

k=l t=l 

Take one of the regular intervals Jq = Jj(fe,m), then there exists to € Jo such that Xfc(O) + T^o = ^m(O). 
Expanding in V(t% — to) we have 

F(x k (0) + Vh) = F{x k {0) + Vto) + F'(x k (0) + Vt )V(tx - t ) + F"(x k (0) + Vt + i9)V 2 (h - to) 2 (37) 

for some < ■& < V(t\ — to). After integration over Jo the linear term in the right-hand side of ((37)) vanishes 
due to symmetry w. r. t. to, constant terms cancel with the corresponding terms in the last sum in (|36[) . and 
as a result we get 

| f (F(x k (0) + Vh) ~ Fix^OWhl <C(F,2)A 3 

J Jo 

The sum of such terms (for given k) over m and over all Jj(fc, m) has the upper bound 

C{F,2)A 2 t (38) 

Taking into account the factor we have the same bound after summation over k. 

Consider now non-regular intervals, For the initial intervals we write the expansion similar to (|3"T|) , taking 
to = 0. Then as above, the constant terms cancel and the linear terms, after integration give bound C(F, 1)A 2 . 
The end intervals Jp, if there is some x m (0) £ Jp, are treated similar and give the same bound. If there are no 
such x m (0) € Jp then we consider the union Jp-x U Jp and do the same procedure. As there are not more than 
2N such intervals, the bound will not depend on t. Namely, this give the bound 

^2NC{F,l)A 2 (39) 

Taking all together, namely Lemma [5) ([3"5)l and (|3"9"|) , we get (|3H|) . 



6 Nonlinear integral equations 

We shall prove the Theorem by taking into account linear and non-linear terms which we skipped in the basic 
linear approximation (|18p . Remind that we assume for simplicity u k (0) — and thus Ci(n) — 0. We consider 
Lj = Lj(k,t) as functions of k and t, given y k (t) and denote 



H{n,t) ==^(y k (t)),L j = L j (n,t) = <Z>(L j (k,t)) 



(we will write down them explicitely below). Applying Fourier transform to (|17j) we get the main system of 
equations in Fourier form 

d 2 H(n,t) A dH(n,t) 4a. q/ . „ n . „. . , ,-, . i ? , n 

^ ' + Yldt + ~M ~ C ° S 2lT N )H{n > <} = M <f>(n, t) + M L ^ *) 

Similarly to (|2"3")l we get the system of integral equations for H(n, t) 

H = K(H) + 7) (40) 
with the non-linear integral operators, acting on H, 

K = J2 K i, K i = K J^ K jAH){n,t) = n(n) / e^W-^LjfatJdti 

j=l 1=1,2 •'° 

and the free term (that we have studied above) 



V(n,t) = Y r l( n ) 

1 = 1.2 J ° 



Define the Banach space Bx = Bj^jv of complex (continuous in t) functions b = b(n, t) on {0, 1, N — 1} x 
[0, T] with the norm 

||6||= sup D n \b(n,t)\ 
n,te[0,T] 

where 

1 

A' 

and for n ^ 



D 



D„ = n A ' 



, MA 

It follows that for any n 



sup \b(n,t)\ < \\bl\D- 1 (41) 

*e[0,T] 



and thus for any subset Q C {0, 1, N — 1} 



]T sup |6(n,t)|< U&H^^ 1 (42) 



Note that 

N-l 



, a 

n— 1 



Y,^<cJ^lnN (43) 



Further on we put for any < j3 < 1 not depending on N 

T = T[fi) = PNC£ g > cfSp^N 3 
We shall prove that Banach fixed point theorem defines the unique solution in Bt- 
Lemma 6 There exists 7 > 0, not depending on N and such that for any sufficiently large N the ball 

D(T, 7 ) = {H €B T :||JT||< 7} 
is invariant with respect to any operator K,Kj, Kj i and for any Hi,H 2 £ Bj- for some q < 1 

\\K{H 1 )-K{H 2 )\\<q\\H 1 -H 2 \\ (44) 

In particular, it follows from Lemma that r/ G D(T, 7). It follows that there exists unique solution and 
we can solve the equation (|40[) by the standard iteration. We will prove Lemma separately for each Kj i with 
sufficiently small q's. 



Bounds for L\ Consider two function yk(t) and Zk{t) of k and their Fourier transforms 

Hi(n,t) = $(y k (t)),H 2 (n,t) = <5>{z k {t)) 

Denote 

Uk(Vk) = F(x k (0) + Vt + y k ) 
Using Taylor expansion we have for some Ok £ [zk, Vk], 6'k e [0> ^fc] 

Myfc) - M^fc) = ufe(zfe + (yfe - z*,)) - uk(zk) = u^\z k ){yk - z k ) + ^u k 2) (z k + k ){y k - z k ) 2 = 

= u ( k\v){Vk - *k) + z k u k 2) (e' k )(y k - z k ) + \u {2 \z k + e k )(y k - z k f (45) 
We shall consider firstly the linear term 

uf{$){y k - z k ) = F^\x k (0) + Vt)(y k - z k ) 

Using the convolution formula 

$(3iff2, n) = ($(5i)*$(52))(n) = 53 ${gi, n - rn)<f>(g 2 , ni) = 53 $(52, n - rci)$(sn, m) 



we have for H = Hi — H 2 



K hl {H) = gn(n) f e*^*-*^ Vx(n - mM)H{mM)dti 

where 

*( n '*) = ^E F(1) ( Xfc (°) + W)exp(2™-) 
fc=i 

Lemma 7 

| X (0,i)l<AC(F,2) (46) 

and /or any < n < 

| X (n,t)| <cn~ x (47) 

Symmetrically for -y < n < N. 
Proof. As for any t 

[ F {1 \x + Vt)dx = 

>/ Si 

we have 

AT 

| X (0,i)l = lx(0,i)- f FW(x + Vt)dx\=\-j-J2FW(x k (0) + Vt)- f F« (x + Vt)dx\ < j-max\F^ (x)\ 

Si k—l 

For n 7^ we can use the following summation-by-parts formula 

JV-l JV-l 

53 h 1 {k)Vh 2 {k) = h^(N)h 2 {N) - /ii(0)/i 2 (0) - ^ + l)V^i(fc) 

fc=0 fc=0 

where fti(iV) = fti(0), /12C/V) = h 2 (0). Put in our case 
Then 

Vh 2 (k) — exp(27rm— ) 



If n < <rN, where for example a = j-, then 

i at nip i ") 
— I h 2 {k + l)Vh 1 {k)\ <max|fca(fc+l)V/ii(fe)| <c 



fc=0 

If crA^ < ?i < then similarly 

1| V* 2 (fc + l)VM*)l<n^M* + l)VM*)l< sup (|- =cC(F,l)n- 1 

N ^ fe l-exp(27rzx) N 

To end the proof we shall prove a general assertion which will be used also in further estimates. 
Lemma 8 If the function x{n,t) satisfies the bounds d^6| ) and d^7| j then under conditions Ii9l 10\) the operator 

K(H)= Y,9n(n) f*e''W t - t ^x{n-ni,t 1 )H(n 1 ,t 1 )dt 1 

1=1,2 "i 

satisfies the bound with q sufficiently small. 
Proof. As in (f2U)) we have by definition 

\\K U {H)\\ < sup D n \gn(n) f e »«W(*-*i> V* jj( n _ ni) tl ) X {n x , t x )dti\ 

n.t Jq 

' u ni 

Firstly, for / = 1 consider, under sup n , the case n = 0, namely 

supAM(O) f e *»(°)(*-*0 Vi?(-m,t 1 )x(ni,ti)^i| < 

o A/f 

<A-V--(£ D- i 1 || J ff||sup|x(n 1 ,t 1 )|+£> o - 1 || J ff||sup| X (0,t 1 )|)< 



< A-^C^IIfflKcJ— +C(F,2)A 2 ^ 

a V ct A \ ai\ 

For the case n / Owe have from (f2~il) . similar to (|27p 



MA 7i V a f—' tl f—' fl 

niT^O "15^0 



The sum of the terms with 71 1 7^ n does not exceed cy In JV, and the term with n\ = n 7^ does not exceed 
—A. Finally we get the bound 

, I MA 1 . . . „ 

c(\/ hxN + -A)g\\H\\ 

V a 71 

If Z = 2, the case n 7^ is similar to above. For the worst possible term with n = 0, 2:2(0) =0we have 
su P C | 5 r 2 (0)|| I dt x ^H{-ni,tx)x(ni,tx)\ < A_1 ^ E sup *)l < 

^A-^t^supixC-m,*)^- 1 !!^! 

The term 



A — t 



lA-^tsuplxl-O^)!^ Il^ll < 2C A , g C(F,2)tA\\H\\ 

A t 

with ni = is small as 

C A , g At < c(3 

The sum of the remaining terms (with m 7^ 0) 

\ A ~ 1 ^X t E ^ V \ X (-ni,t)\D-}\\H\\ < 2C A , g J2^]f^\\H\ 

ni^O ni 1 

is small by the right inequality of ©. 



The second (nonlinear) term 

z ku\ '(0 k )(y k - z k ) 

is intuitively simpler because of the additional small factor z k . Using the Fourier transform \i *H(n) is treated 
similarly to the first term, where instead of x( n ) t) one should take 



Xi(n,t) = ZkU< k ] ( e 'k) ex P( 2 ™™^) =X2*H 2 

k=l 

where 



Xa (n,t) = I^4 2 )(^)exp(27rm^) = i £ ^ (2) (^(0) + W + <£) exp(27rm|) 
fc=i fe=i 

We will need only the obvious bound 

|X2 (n,t) | <C(F, 2) (48) 
but we can prove more. In fact, we have by (|4"2")) and (j4Tfl) 



Bup| VJb (i)| < V sup |fli(n,t)| < \\HxW Vd; 1 < HflilKA + cJ— InJV) < c(l + p)A (49) 



as © gives 



— lniV<cpA (50) 



The same holds for \z k \ and then also 

\e k \,\e' k \<cA 

This shows that the estimates of Lemma [7] hold also for X2( n j^)j with the same proof. 
Then by (@2J and ((5DJ) 

IXi (n, «)| = |X2 * ffal (», *) < C(F, 2) £ |# 2 (n, t) | < C(F, 2) | |ff 2 1 1 ]T D7 n l 

n 

Then using lemma [HJ we get the result. The third term in (|45[) is treated similarly. 

Bounds for L2 We have 

L 2 = 2aA- 3 V-[<^,Vy fc ] = 2aA- 3 ^V-Vy fc + 2aA- 3 (V-^)(S*- 1 Vy fe ) 

L a (n, i) = 2aA- 3 Ci, 2 * $(V" Vy fc ) + 2aA- 3 Ci,i * ^S^V^), 

where 

1 w 

Ci >2 (n) = $ (-5 1 ) = ^ E ^ ex P( 2 ™™^)> Ci,i(n) = ^(V-5 1 ) 
fc=i 

Lemma 9 For any n 

ICmWI < ca^gA 2 (51) 

|Ci,2(n)| < ca^gA (52) 
This easily follows from Lemma [TJ series ([5]) and bound (f55"j) . 
Lemma 10 For g = 1,2 f/ie operators 

K n<q {H)=n(n) f e*'W( t - tl ^(exp(-27ri^^)-l)«fr(n-ni > ii)Ci, g (ni,ti)d«i 

•'O ni 

satisfy the bounds 

WK^H^Ka^gAWHW 



Proof. We shall demonstrate the (straightforward) calculation for q = 2. The case q = 1 is quite similar and 
even easier (because of A 2 factor). 

For 1 = 1, as in (|2U|) we have by definition and by (ptTl) 

11^2,1.2(^)11 <sup£>„|n(n) /^^W^) V(exp(-2^^— ^)-l) 2 ff(n-n l! t 1 )C 1 , 2 (n 1 ^ 1 )di 1 | 

Consider, under sup„, first the case n = 0, namely 

Bupi?o|ri(0) /"VOX'-'*) V(exp(27ri^)-l) 2 ff(-T ll ,i 1 )Ci.2(ni,*i)d«i| < 



<A- 1 |^(E^i 1 ll ff ll^^( B i.*0<4T l|ff||a "W^ 1,lJV 

"i#0 1 



After multiplying on aA 3 we get 



1M \iTTU -1 / MA M """ /M 

ll l|I||fl 5 V~ 

For the case n/0we have, similar to (|2~T|) 



a ~ 3 aT i|h||q_15 V ~ lniV - cA_2c ^T ||i;r|| VS lnAr - cp||ff| 



< aA- 3 ||iJ||AcC* Q , g A < csA -1 ||i2"|| 

Similarly for I = 2 

||#W*0ll <su Pj D„|r 2 (n) f e ^(n)(t-ti) V(exp(-2™^— !*) - f) 2 ff(n - m, ti)Ci. 2 (rai, *i)*il 
„,t Jo ^ AT 

The case n 7^ is similar to above. For the worst possible term with n = 0, £2(0) = we have, noting that the 
term with n% = is zero, 

sup A) MO) I |i V(exp(2™^) - lfH{- ni MXi,%{nuh)\ < A" 1 ^ Y\ sup |Ci,a(-»i, t)\D^\\H\\ < 
m N A ut#o ( 

< A _x 4*ca -1 ffAcpA||fr|| < cpita _1 gA||ff|| (53) 

as 



^C 1 < c\I^—\d.N < cA\ ^r-hxN < cpA 
*—f \ a V aA 

m#0 



Finally 



aA- 3 cp-to- 1 5 A||ff|| < cptA- 2 C A J\H\\ < cp(3\\H\ 



1 

A 

The last bound follows from (|10[) . The lemma is proved. 

Bounds for L3 We shall consider the m-th term L^ m of L3 in the series (JT6 



A m+2 ~ u a m+2 \ v A v yK\")) ' "V" v yK^jj \ v ; . m+2 

The convergence of the series in m will follow from the obtained bound for L^. m . Using the Leibnitz formula 
(fTl|) we can rewrite the first term as follows 

1 1 m 

A fe A /c , = 1 



Its Fourier transform will be 

<i> («T^+2(-V-)(V + ?A W) m )=^Ci(r i -n 1 ,i)<I'(V-V + t/ fe (t))(n) 



where 

1 m 

Ci(M) = *(-«^ pi E( V+ ^(*)) ? '( 5 " 1V+ ^(*)) ro ^'" 1 ) 



AS 

k 2=1 



Using Lemma |U we have for m > 2 



|Ci(n.*)l < aA- m - 2 m(cC Ag \ — In TV)™- 1 < aA-™- 2 ™!^ 4 )" 1 - 1 < aAmfcp)™- 1 

a 



Also by 

|v-v+tfr(t)| < 2 l#(Ml < E^ 1 !!^!! ^ c <° A H ff H 

Then we have the result similarly to (f51?l) . The second term is treated similarly. It is interesting to note that 
nonlinear terms demand less restrictive bound than (|10[) . 

7 Comments 

There are many problems left. 

1. Most irritating and interesting is however only one: to include the ionic lattice to the model of strongly 
interacting electrons. May be a satisfactory picture can be obtained only on the quantum level. However 
on the quantum level it is not clear even how to write down the Schroedinger equation because the external 
field F(x) is not potential on the circle. 

2. With our methods we could not prove stability for any time t € (0, oo) because of the zero mode problem, 
that is existence of zero root for n — 0. Additional linear term proportional to yk in the basic equations 
(fT3|) could easily solve this problem but I could not obtain this term as a result of realistic interaction 
with the ionic lattice. 

3. Our assumption concerning smallness of g is too restrictive at least in two points. Firstly, if a~ x g is 
bounded then the space scale TV -2 controls the effective forces acting on the electrons. If a is smaller 
than g, then this scale will be in-between A -1 and N~ 2 , but when it becomes comparable with the scale 
A -1 , then the macroscopic homogeneity will be lost. In particular, the macro- velocity V may depend on 
the distance from the support of the external force. 

4. Secondly, the worst perturbation term is the linear term Li . Possibly, more refined techniques allow better 
estimates. 
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